function diff = q_e(m,psi,param)
% q in the replacement region

    diff = q_0_e(m,psi,param) + q_1_e(m,psi,param);

end
